Ultracold atoms in optical lattices generated by quantized light fields 
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We study an ultracold gas of neutral atoms subject to the periodic optical potential generated 
by a high-Q cavity mode. In the limit of very low temperatures, cavity field and atomic dynamics 
require a quantum description. Starting from a cavity QED single atom Hamiltonian we use different 
routes to derive approximative multiparticle Hamiltonians in Bose-Hubbard form with rescaled or 
even dynamical parameters. In the limit of large enough cavity damping the different models agree. 
Compared to free space optical lattices, quantum uncertainties of the potential and the possibility 
of atom-field entanglement lead to modified phase transition characteristics, the appearance of new 
phases or even quantum superpositions of different phases. Using a corresponding effective master 
equation, which can be numerically solved for few particles, we can study time evolution including 
dissipation. As an example we exhibit the microscopic processes behind the transition dynamics 
from a Mott insulator like state to a self-ordered superradiant state of the atoms, which appears as 
steady state for transverse atomic pumping. 



I. INTRODUCTION 

Laser light, far red detuned from an atomic resonance, 
is nowadays a standard tool in experimental quantum op- 
tics to create tunable optical potentials [ij which can be 
loaded with ultracold atoms to provide for a plethora of 
possibilities to study quantum properties of many-body 
strongly correlated systems Q ■ The high level of micro- 
scopic understanding and extensive control of the light 
fields and atoms allow to implement genuine models like 
e.g. the Bose-Hubbard (BH) model [3|Jj|. Initially origi- 
nating from condensed matter physics [5| it has been used 
to study the Mott insulator to superfluid phase transi- 
tion [|| in detail and in real time. Adjusting several of the 
lattice parameters as the intensity and the configuration 
of the lattice lasers provides a versatile toolbox of tech- 
niques to control the dynamics of the atoms in the lat- 
tice [7[. Moreover, the collisional properties of the certain 
types of atoms can be tailored by means of magnetic @ 
or optical 0] Feshbach resonances. Using extra confine- 
ment it was even possible to observe the Mott insulator to 
superfluid transition in ID [l(| [H[ and 2D [12] , followed 
by other spectacular demonstrations of condensed matter 
physics phenomena as the realization of a Tonks gas in 
ID [ID, [3] and the Berezinskii-Kosterlitz-Thouless phase 
transition in 2D fl5| . Theoretically many more proposals 
to apply these methods to spin systems and investigate 
further fascinating properties of strongly correlated sys- 
tems were put forward (see 16] for a review). 

In all of these approaches, the light fields were approx- 
imated by classical, externally prescribed fields indepen- 
dent of the atoms. This requires intense light, far detuned 
from any atomic transition. Of course this assumption 
holds no longer if the light, which generates the optical 
lattice, is enhanced by an optical resonator. In this case 
- given a sufficient atom number N and atom-field cou- 
pling g - the field itself becomes a dynamical quantity (l7j 
depending on the atomic distribution. As all atoms are 
coupled to the same field modes, this immediately intro- 
duces substantial long range interactions, which cannot 



be ignored as in free space. In specially designed cases 
this force induces coherently driven atoms to self-organize 
in regular patterns as predicted in Ref. [H, and sub- 
sequently experimentally verified (2C| . 

In addition, in a high-Q optical resonator relatively low 
photon numbers are sufficient to provide strong forces. 
This was demonstrated by trapping an atom in the field 
of just a single photon [2l|, [22|. Hence the inevitable 
photon number fluctuations induced by cavity damping 
generate force fluctuations on the atoms causing diffu- 
sion. At the same time as cavity photon loss constitutes 
a dissipation channel, it can also carry out energy and en- 
tropy of the system. This opens possibilities for cooling of 
atomic motion [23|, HJ, [2f| [26[ , as demonstrated by beau- 
tiful experiments in the group of Rempe [U [28| . Since 
this cooling mechanism does not require the existence of 
closed optical cycles it could even be used for qubits [29[ 
or to damp quantum oscillations or pha se fluctuations of 
a BEC coupled to a cavity field [30l. l3l| . 

For low photon numbers the quantum properties of 
the light field get important as well and the atoms are 
now moving in different quantized potentials determined 
by the cavity photon number. Quantum mechanics of 
course allows for superpositions of photon numbers in- 
voking superpositions of different optical potentials for 
the atoms. First simplified models to describe this new 
physics were recently proposed by us [32j and in paral- 
lel by other authors [33]. As the intracavity field itself 
depends on the atomic state (phase), different atomic 
quantum states are correlated with different states of the 
lattice field with differing photon number distributions. 
In this way quantum mechanics allows for the creation 
of very exotic atom-field states, like a superposition of a 
Mott-insulator and superfluid phase, each thereof corre- 
lated with a different photon number. Some quite exotic 
looking phase diagrams for this system were already dis- 
cussed in Ref. [33]. Without resorting to the full com- 
plex dynamics of the system, the quantum correlations 
between the field and the atomic wavefunctions open 
the possibility of non-destructively probing the atomic 
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state by weak scattering of coherent light into the cavity 
mode 34] and carefully analyzing its properties (3f|. 

It is quite astonishing, that experimental progress in 
the recent years has made such systems experimentally 
accessible and at present already several experimental 
groups succeeded in loading a BEC into a high-Q optical 
cavity [H, [13, El, El, Ej. A reliable analysis of these 
experiments has made more thorough theoretical studies 
of such systems mandatory. 

In this work we concentrate on the study of an ultra- 
cold gas in optical lattices including the quantum na- 
ture of the lattice potential generated from a cavity field. 
This extends and substantiates previous studies and pre- 
dictions on such a system by us [HJ as well as other 
authors [Hj]. Here we limit ourselves to the case of a 
high-Q cavity which strongly enhances a field sufficiently 
red detuned from any atomic transition to induce an op- 
tical potential without significant spontaneous emission. 
In particular we address two different geometric setups, 
where either the cavity mode is directly driven through 
one mirror, or the atoms are coherently excited by a 
transverse laser and scatter light into the cavity mode. 
The cavity potential can also be additionally enhanced by 
some extra conservative potential applied at a different 
frequency [4l|, [42[ ■ These two generic cases leads to quite 
different physical behavior and allow to discuss several 
important aspects of the underlying physics. 

This paper is organized as follows. Sec. HI] is devoted 
to a systematic presentation of our model and various 
simplifying approximations as adiabatic elimination of 
the excited states of the atoms and subsequent formula- 
tion of an effective multi-particle Hamiltonian in second 
quantized form. In section HTT1 we specialize on the sim- 
plest generic case of a coherently driven cavity and ap- 
proximate the corresponding Hamiltonian by adiabatic 
elimination of the cavity field. We investigate the prop- 
erties thereof, corresponding to the influence of the cav- 
ity on the Mott-insulator to superfluid quantum phase- 
transition and identify the regime of validity for the elim- 
ination of the cavity field. Finally, we compare these re- 
sults with the dynamics of the full master equation. In 
Sec. HVI we study the more complex case of atoms coher- 
ently driven by a laser field transversal to the cavity axis, 
where it is much harder to find valid analytical simplifi- 
cations and one has to resort to numerical studies of few 
particle dynamics. Finally, we conclude in Sec. fVl 



II. MODEL 

We start with N two-level atoms with mass m and 
transition frequency uj eg strongly interacting with a single 
standing wave cavity mode of frequency u> c . We also 
consider coherent driving of the atoms at frequency lo p 
and with maximal coupling strength ho and of the cavity 
with amplitude i] (see Fig. [T]). Note that in the specific 
examples later we will consider only one pump laser beam 
at a time. 
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FIG. 1: (color online). Scheme of atoms inside an optical 
cavity, driven by two external pumping lasers. An additional 
conservative lattice potential, independent of the intracavity 
field, is realized by a far off- resonant dipole trap (FORT). 



Using the rotating-wave and electric-dipole approxima- 
tion, we can describe a single atom of this system by the 
Jaynes-Cummings Hamiltonian [43| 
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Explicitly the different Hamiltonians for the atoms, the 
field mode and the interaction read: 

-2 
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2m 

- ihh{*) (a+e- 1 ^ ~ a-e Wpt ) , (2a) 
h£ ] = hw c a)a - ihq (ae 1 ^ - aV^*) , (2b) 



H 



(i) 

Int 



-ihg(pc) [a a — a a 



(2c) 



Here /i(x) denotes the mode-function of the transverse 
pump field, g(x) denotes the cavity mode function and 
the field operator a describes the annihilation of a cavity 
photon with frequency oj c . V e (x) and V s (x) are external 
trapping potentials for the atom in the excited and the 
ground state, respectively. In order to change to slowly 
varying variables we apply a unitary transformation with 
operator U(t) = exp[iu> p t (a + a~ + a^aj], such that we 
end up with the following single-particle Hamiltonian, 
using the same symbols for the transformed quantities: 
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(3a) 
(3b) 
(3c) 

denotes the detunings 
of the cavity and the atomic transition frequency from the 
pumping field frequency. In order to describe the situa- 
tion for N atoms, we use the single-particle Hamiltonian 
of Eq. ([T]) and in second quantization formalism 44j , 
i.e., 



where A r 



-U) c 



H = H A + H R + Ha-r + H A - P + H A _, 



(4) 
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The terms in this expression correspond to the single 
particle terms in @ and ([3]). Hence, Ha and H R model 
the free evolution of the atomic and the field variables, 
respectively. They read as: 



H A = I d 3 x 



2m 



V 2 + K,(x) * 9 (x 



2m 



V 2 -fiA Q + y e (x) * e (x) 



(5) 



where x f , g (x) and ^(x) denotes the atomic field opera- 
tors for annihilating an atom at position x in the ground 
state and the excited state, respectively. They obey the 
usual bosonic commutation relations 



^(x^UxOl = 6*{x-x?)6 f , f , 



*,(x), = *t(x),*Ux^ 



0, 



(6a) 
(6b) 



for /, /' G {e,<?}. The field operator remains unchanged, 
i.e., Hr = —hA c a^a— ihq (a — e^). The two-body inter- 
action is modeled by a short-range pseudopotential [45| . 
characterized by the s-wave scattering length a s , leading 
to a Hamiltonian 

Ha-a = || d 3 x*t(x)*t( x )vl/ 9 (x)vl/ 9 (x), (7) 

where U = Ana s h 2 jra. The coupling of the cavity field 
with the atoms inside the cavity is given by 



H A -R = -ih J d 3 x*J(x)g(x)a t $ e (x) + h.c. 



(8) 



whereas the interaction with the laser beam, which co- 
herently drives the atoms, reads 



H A -i 



ih J d 3 x*t(x)/i(x)* e (x) +h.c. 



(9) 



Let us now calculate the Heisenberg equations for the 
various field operators, starting with the operator for the 
excited state, i.e., 



dt 



K(x) 



2m 



*e(x) 



[. 9 (x)a + Mx)]* s (x). (10) 



The first term corresponds to the free evolution of the 
atomic state, whereas the second term describes the ab- 
sorption of a cavity photon or a pump photon together 
with the annihilation of a ground state atom. Similarly, 
the equation for the ground state operator reads: 



a* g (x) 

dt 



2m 



-V 2 



ft *J(x)* fl (x) 



+ [g(x)a t +/i(x)] * e (x). (11) 



Finally, the Heisenberg equation for the cavity field op- 
erator is given by: 

^ = iA c a + r, + J d 3 x.g(x)* 3 (x)* e (x). (12) 

Again, the first term corresponds to the free field evolu- 
tion, whereas the last two terms are driving terms of the 
cavity field. 

As we want to treat temperatures close to T — we 
have to avoid heating and ensure weak atomic excita- 
tion, where there is only negligible spontaneous emission. 
In this limit we can adiabatically eliminate the excited 
states from the dynamics of our system. This requires 
large atom-pump detunings A a , where we also can ne- 
glect the kinetic energy term and the trapping potential 
in (fT0|) compared to A a . Necessarily, we assume that the 
field operators 'l'g(x) and a vary on a much slower time 
scale than the 1/A a terms, such that we obtain: 

* e (x,t) =--J-[/i(x)+ fl (x)o(t)]* 9 (x,t). (13) 

Inserting this expression for 'l'e(x) into (|ll|l and (TT2"|) 
leads then to: 



dt 



V g (x) /* 2 (x) 



2 (x) 



2m 



Mx)ff(x) 
A a 



da 
St 



(a + . 



A a 



A a 



n 
u 



a/ a 



j.*t(x)* 9 (x) 



* fl (x), (14) 



d 3 x 5 2 (x)M/t(x)* s (x) 



d 3 x/ l (x)*t( x )^( x)+7? . (15) 



To discuss the underlying physics in a tractable form, 
the trick is now to find an effective Hamiltonian H c g 
which leads to the same dynamics as given by Eq. (fT4]) 
and (|15p . Thus this Hamiltonian has to obey: 

lh ^dt^ = -Heff] and ih^ = [a,H eS ] . (16) 

From this we can easily read off a possible effective Hamil- 
tonian of the form: 



H e s ■ 

h 



d 3 x*j.(x) 



2m 



v 2 + y 9 (x) 



[/i 2 (x) + 3 2 (x)at a + ft(x) 5 (x) (a + a 1 )] j * fl (x) 

+ || d 3 x*t( x )vl,t(x)v|j 5 (x)* ff (x) 

— ihr) (a — a^) — hA. c a) a. (17) 

The corresponding single particle Hamiltonian, which 
leads to this second quantized Hamiltonian is [46] ]: 

H$ = + K 9 (x) + A [ft» (x) + ff 2 (x)a t a 
+/i(x)s(x) (a + a f )] - i^?7 (a - a f ) - ^A c a f a. (18) 
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This simplified effective atom-field Hamiltonian will be 
the basis of our further considerations. It is, however, 
still much too complex for a general solution and we will 
have to make further simplifying assumptions. Hence at 
this point we will restrict ourselves to ID motion along 
the cavity axis. In an experimental setup this could be 
actually realized by a deep radial trapping potential, but 
we think that at least qualitatively the model should also 
capture the essential physics if some transverse motion 
of the particles was allowed. As one consequence this 
assumption requires a rescaling of the effective two-body 
interaction strength [47| , which enters as a free parameter 
in our model anyway. 

Mathematically we thus end up with a one-dimensional 
optical lattice, wich is partly generated by the resonator 
field and superimposed onto a prescribed extra trapping 
potential V^(x) = V g (x). The mode function of the cav- 
ity along the axis is approximated by <?(x) = g(x) = 
go cos(kx) and the transverse laser beam forms a broad 
standing wave h(x) — hocos(k p y), which in our one- 
dimensional considerations (y = 0) is just a constant 
term that we can eventually omit in (|17p . 

As we consider external pumping of atoms and mode, 
we essentially treat an open system and we have to deal 
with dissipation as well. Such dissipation processes are 
modeled by Liouvillean terms C appearing in the master 
equation for the atom-field density operator, i.e., 



9= ~z [H e e, q] 
in 



Cg. 



(19) 



As mentioned above, we assume large atom-pump de- 
tuning A a , suppressing spontaneous emission to a large 
extend. However, we still have to deal with the cavity 
loss k, which will thus be the dominant dissipation pro- 
cess. Hence the corresponding Liouvillean using a stan- 
dard quantum optics approach [50j reads: 



Cg = k (2aga) — a) ag — ga'aj 



(20) 



Equivalently in the corresponding Heisenberg equation 
for the field operator, cavity loss leads to damping terms 
and fluctuations, so that it then reads: 

A c - j^- J dx^l(x)cos 2 (kx)^ g {x) - k\ a 
_ t 9^o f dx ^ g ^ C0S ^ kx ^ g ^ +ri + T . n _ (21) 

Since we will be mainly interested in normally ordered 
quantities and assume vacuum (T=0) outside the cav- 
ity, the input noise operators Ti n will not enter in the 
dynamics, such that we will omit them later. 

Let us now proceed and transform the Hamiltonian 
into a more commonly known form. Following standard 
procedures, one constructs maximally localized eigen- 
functions at each site and expands the atomic field oper- 
ator ^> g (x) in terms of single atom Wannier functions [49[ 



where b n ,k corresponds to the annihilation of a particle 
in the n— th energy band at site k. Since we assume the 
involved energies to be much smaller than the excitation 
energies to the second band, we are able to keep only the 
lowest vibrational state in the Wannier expansion, i.e., 
^g( x ) = J2kbkw(x - x k ), where w(x) = w (x). This 
yields to the following Hamiltonian: 



H 



^g( x ) = X! X! b ri,kW n {x - x k ), 



(22) 



Ekiblh + (fitWa + Vol) Jklb l bl 

k,l k,l 

+ KrjeS (a + a)) ^ Jklb\h - ifvq (a - a)) 

k,l 

+ \ J2 U t]kl b\b)b k h-h^a, (23) 

i : j.k.l 

where the addendum eff of the Hamiltonian is omitted. 
Here we introduced an important characteristic param- 
eter of atomic cavity QED, namely the refractive in- 
dex Uq of a single atom at an antinode, which is given 
by Uq = <7o/A . It gives the frequency shift of the 
cavity mode induced by a single atom at an antinode 
and also corresponds to the optical lattice depth for an 
atom per cavity photon [17] . Similarly, the parameter 
Vctt — 9o^o/A a describes the position dependent effec- 
tive pump strength of the cavity mode induced by the 
scattered light from a single atom at an antinode. 

Note that the Wannier state expansion Eq . |2"21 depends 
on the potential depth. Thus the Wannier functions and 
the corresponding matrix elements depend on the cav- 
ity field and thus in principle are dynamic quantities. 
However, they keep the same functional form with a few 
changing parameters, which have to be determined con- 
sistently. This is of course consequently also true for the 
various coupling parameters in the Hamiltonian. The 
above model thus can only be valid as long as the single 
band approximation stays valid during the system dy- 
namics and the parameters dont change to rapidly. In the 
special but rather interesting case, where the atoms are 
trapped solely by the cavity field [U [22j this is not valid 
for very low photon numbers. Here a single photon num- 
ber jump will induce excitation to higher bands, which in- 
duces nonlinear dynamics beyond the single band model. 

In practise this problem can be circumvented by adding 
an additional external trapping potential V g (x) to the 
model, which guarantees a minimum potential depth 
even in the case of zero cavity photons. Experimen- 
tally this is feasible, for instance, with a far detuned, 
off- resonant dipole trap (FORT) [Fl|, i.e., V g (x) = 
V c iCOS 2 (kpx), where kp denotes the wave number of the 
FORT field. In the experimental realization, the fre- 
quency of the corresponding laser field ujp is only very 
few free spectral ran ges separated from the main cav- 
ity frequency uj c [271 . 152L [53|. Hence, in the vicinity of 
the cavity center, the coincidence of the FORT field and 
the cavity field is very good, and we can replace in good 
agreement cos 2 (fc_F^) with cos 2 (fcir). 

Let us remark here that by including this extra poten- 
tial, we can keep our model and allow for further analyti- 
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cal analysis of the dynamics, but we also have thrown out 
a great deal of interesting physics already. Actually, for 
very few atoms one still can solve the full Hamiltonian 
without the restriction to the lowest bands by quantum 
Monte Carlo wavefunction simulations. Some early re- 
sults of such simulations can be found in Ref. 0, 
However, this is not the subject of this work and we will 
proceed here with the effective lattice model under the 
assumption of a deep enough extra potential or strong 
enough cavity fields. 

Note that in (j2li|) . in contrast to the case of the Bose- 
Hubbard model in a classical optical lattice, where the 
matrix elements of the potential and kinetic energy can 
be merged, here two separate parts exist due to the pres- 
ence of the cavity held operators in the Hamiltonian. Ex- 
plicitly they read as: 



E u = / dxw(x - x fc ) ( -7^;V 2 ) w(x ~ xi), (24a) 



JfcZ 



•hi 



dx w(x — Xk) cos 2 {kx)w{x — xi), (24b) 
dxw(x — Xk)cos(kx)w(x — xi). (24c) 



The on-site elements Jkk and Ekk are independent of the 
lattice site k, whereas Jki changes sign periodically, i.e., 
Jkk = —Jk+i,k+i due to the cos, which has twice the peri- 
odicity of the lattice. This also accounts for Jk,k+i = 0. 
Note that the existence of this term implies that two 
adjacent wells acquire different depths forcing us to re- 
assure that for the case of the directly pumped atom 
rj e ff(a + a^) cos(kx) is only a small perturbation of the 
lattice. As the next-nearest elements are typically two 
orders of magnitude smaller than the nearest-neighbor 
term [3[ they can safely be neglected (tight-binding ap- 
proximation). Hence we label the site- independent on- 
site matrix elements with Eq , Jo and Jo , whereas E and 
J are the site-to-site hopping elements. Furthermore, in 
the case of the nonlinear interaction matrix elements, 



Uijki = giD J dxw(x - Xi)w{x - Xj)w(x - x k )w(x ~ xi) 

(25) 

we can omit the off-site terms since they are also typ- 
ically two orders of magnitude smaller than the on-site 
interaction matrix elements. Note that gio is the one- 
dimensional on-site interaction strength, originating from 
an adjustment of the scattering length a s , due to the 
transversal trapping [47|. As a central result of our 
studies we therefore obtain a generalized Bose-Hubbard 
Hamiltonian: 



H 



= E N + EB+ (hU a)a + V c i) (j a N + JB 
+ hrj cS (a + a f ) J y^(-l) fc+1 n fc - HA c a^ a 



U 



ihr) (a - a f ) + — h k (h k - 1 



(26) 



where the nonlinear on-site interaction is characterized 
by U = gin J dx |u>(ai)| 4 . In addition, we introduced the 
number operator N — ^ fe hk = J\ ^l^fe and the jump 
operator B = ^ fe ffyt+i&fc + h.c^j . Note that for strong 
classical intracavity helds and no transverse pump we 
recover the standard Bose-Hubbard Hamiltonian. 

Finally, let us remark that we now can also rewrite the 
held Heisenberg Eq. (f2"Tj) in the above terms, which gives: 



(j 



N ■ 



JBj^ - k| a + 77 
-^ off Jo^(-l) fe+1 ^. (27) 



Here we clearly see that besides the number operator 
N for the atoms also the coherence properties via the 
operator B and statistics via hk play a decisive role in 
the held dynamics. As this held acts back on the atomic 
motion, interesting and complex coupled dynamics can 
be expected from this model, which was partly already 
discussed in [32|, [33| and will be elucidated more in the 
remainder of this work. 



III. CAVITY PUMP 

Let us now turn to the conceptually simplest case and 
restrict the pumping only to the cavity, where only a 
single mode is coherently excited (cavity pumping). This 
mode will generate an optical potential in addition to the 
prescribed external potential. For large enough photon 
numbers the external potential can even be omitted and 
the particles are trapped solely by the cavity held. As es- 
sential ingredient in the dynamics, the identical coupling 
of all atoms to this same held mode induces a long-range 
interaction between the atoms independent of their po- 
sitions. Setting 7/ ff = 0, the Hamiltonian (|26]) is reduces 
to: 



H = E N + EB+ (hUoa^a + V d ) (j N + JB 



hAcC^a — ihrj (a — a^) + —C '. 



(28) 



Here wc introduced C = h k (hk — 1) for the operator 
of the two-body on-site interaction. Still we see that the 
corresponding Heisenberg equation for the cavity field: 



|i A c - U ( J N + JB} -k\ 



(29) 



depends on atom number and coherence. For very weak 
fields this yields an atom statistics dependent cavity 
transmission spectrum, which was studied in some de- 
tail in Ref. [35j. Here we go one step further and study 
the dynamical back action of the field onto atomic motion 
and field mediated atom-atom interaction, which appear 
at higher photon number. As the model is still rather 
complex we need some further approximations at this 
point in order to catch some qualitative insight. 
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A. Field-eliminated Hamiltonian 

Although the influence of the cavity field on the atoms 
is equal on all particles, their common interaction gen- 
erates a dynamics much more complex than for a Bose- 
Hubbard model with prescribed external potential. This 
is more analogous to real solid state physics where the 
state of the electrons also acts back on the potentials. To 
exhibit the underlying physics, we will now derive an ap- 
proximate Hamiltonian, which solely depends on particle 
variables by adiabatically eliminating the field |28|) . This 
should be valid when the damping rate K of the cavity 
generates a faster time scale than the external atomic de- 
grees of freedom. Actually as tunneling is mostly a very 
slow process (much slower than the recoil frequency) , this 
will be almost always the case in practical experimental 
setups. To this end, we simply equate (f2U)) to zero and 
obtain formally a = r)/{n — i[A c — Uq(JoN + JB)]}. In 
the following we constrain ourselves to the case of a fixed 
number of atoms, i.e., N — Nl. The very small tun- 
neling matrix element J can be used as an expansion 
parameter, leading to: 



n 



iA' 



U Q J 



iA' 



■B 



(UoJ) 2 
(k - iA' c 



:B< 



(30) 



where we introduced a shifted detuning A' c — A c — 
U Q J N. 

In order to obtain an effective Hamiltonian, where the 
cavity degrees of freedom are eliminated, we replace the 
field terms in |28|) . by the steady state expressions (|30"|) . 
as well as in the Liouville super operator (f20|) . Note, 
that this is more appropriate than the naive approach of 
a replacement just in the Hamiltonian, as has been done 
in our former work [32[ . If we consider terms up to order 
cx J 2 , the exchange in the Hamiltonian yields: 



ffad = {E+JV d )B+jC 



A'' 



?B 



(31) 



Altogether, we end up with a Hamiltonian, where the 
cavity field has been eliminated: 



H. 



ad 



= (E + JV cl )B 
W^Jr, 2 1 



U 



c 



B 



2 

u JK 



3A'/ 



B 2 



(33) 



The loss rate of the cavity is described by the remaining 
dissipative part of (|3"2l : 



C&dQ = 



K U 2 J 2 V 2 
(n 2 +A' 2 )' 



2BgB - B q — qB' 



(34) 



Note, that the above adiabatic elimination procedure 
is not completely unambiguous due to ordering freedom. 
Nevertheless it should give a qualitatively correct first in- 
sight. An alternative way of deriving an effective Hamil- 
tonian, depending solely on particle observable is similar 
to ([in]) and (fl~7|) . This amounts to a replacement of the 
field variables with (pjU)) in the Heisenberg equation for 
the external atomic degrees of freedom, which read as 
follows: 

hk = \ i( E + JVcl + hu oJa f a) (bk-i + bk+i) - Uh k b k ] . 

A naive replacement of the field operator a and its ad- 
joint a' by (|5D)l in the above expression leads to an equa- 
tion for 6^, which cannot be generated from an effective 
adiabatic Hamiltonian in the form bk = — i/h[bk, -ff a d]- 
Hence, before substituting the adiabatic field operators, 
we have to symmetrize the expression containing the field 
term in (|35[) in the form 



b k = --[(E + JV cl ) (6 fc _! + 6 fe+ i) - Un k b k ] 

- [a f a (bk-i + bk+i) + (&fe-i + b k +i) a) a] 



(36) 



Next, by applying the same procedure to the Liouville 
equation - again up to terms cx J 2 - we obtain an adia- 
batic Liouville operator: 



£ad£> : 



2U Jn 2 V 2 

(« 2 + A^ 2 ) 2 
kU 2 J 2 V 2 
' ( K 2 + A> 2 ) 2 



B 



2A' C U J 



(2BgB 



- B 2 g- qB 2 



(32) 



The Lindblad terms in the second line are real, cor- 
responding to dissipation, whereas the first, imaginary 
term corresponds to a unitary time evolution and has 
therefore to be added to the adiabatic Hamiltonian, i.e., 



H, 



ad 



H, 



ad 



2hU Jn 2 V 2 
U 2 + A' 2 ) 



B 



2A' c U a J 
K 2 + A' 2 



B' 



This form enables us to describe the dynamics of bk 
by a Heisenberg equation with an effective Hamiltonian, 
which up to second order in J reads: 



tfad = {E+ JV cl )B + 

hUoJr/ 2 



B 



(37) 



The terms in the second line stem from the field terms 
in (|36|) . Although this Hamiltonian looks a bit different 
from the first version derived before (|33[) . their properties 
are - within their regime of validity - in very good agree- 
ment as long as hopping is slow compared to damping. 

To exhibit the physical content of this Hamiltonian 
one can look at its eigenstates. As first step we calcu- 
late the Mott insulator state [see Eq. ||4T)|) ] fraction of the 
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lowest energy state \ip) of these two Hamiltonians, i.e., 
Pmi = I (V'IMI) | 2 (see also Fig. 2]), as a function of the on- 
site interaction energy for different values of A' c . This will 
indicate changes of position and behavior of the Mott in- 
sulator superfluid transition (see Fig. [4]). To compare the 
two approximate Hamiltonians in Fig. [2j we plot the dif- 
ference of the Mott insulator fraction of the ground state 
of (|33p and (|3"T)) , as well as the difference of the steady 
state photon number. Obviously the two Hamiltonians, 
converge in the limit of large cavity decay k. This can 
also be seen in Fig. ® where the dashed-dotted line de- 
picts the case of a smaller A' c (which is equivalent to an 
enlarged k), showing a strongly enhanced coincidence. 




0.05 0.5 1 1.5 

g 1D /(E R d) 



B. Field-eliminated density operator 

Let us now use a further and somehow more systematic 
alternative approach to eliminate the cavity field dynam- 
ics from the system evolution directly from the Liouville 
equation by following a method proposed by Wiseman 
and Milburn [54j , which is valid for large n and low pho- 
ton numbers. In this case we have 



(Hat) 




hUo(a^a) 


K 




K 



= e <C 1, 



(38) 



where H at is the atomic part of (|28|) . i- e -, H a t = 
(E + V ci J) B + UC/2. Again the total atom number N 
is supposed to be constant. This allows to expand the 
density operator in powers of e, corresponding to states 
with increasing photon number: 

Q = Qo® |0)„<0| + (ei® |l>«(0| + ft.c.) (39) 

+ 92 ® |l)a(l| + (92 ® |2)a(0| + h.C.) + O ( £ 3 ) . 

Here Qi are density operators for the particle variables, 
corresponding to the order i of magnitude in the expan- 
sion parameter e. We substitute this expression into the 
Liouville equation (fT9j) with the Hamiltonian from ([28]) . 
which leads to the following set of equations: 



go 

Qi 

i>2 



■4 [H at , Qo] -v(qi + Qi) + 2k £>2 



(40a) 



ih 



[H at , gi] - i] (V2g' 2 + q 2 - go) - K£i(40b) 



+ i 



A c -U { J N + JB 



I 

ih 
iU 



qi + kO (e 4 ) 
[Hat , Qi] + V ( Qi + g\ ) - 2k£» 2 (40c) 



J a N + JB, g 2 



g-2 = -k[H a uQ' 2 ]+V2r)Qi-2KQ 2 
in 



(40d) 



2i 



A c -U J N + JB 



g' 2 + kO (e 4 ) 



Now we adiabatically eliminate the off-diagonal elements 
Q\ and £>2- Setting their derivations in (|40b|) and (|40d[) 




0.05 0.5 



1 1.5 

9id/ {E R d) 
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FIG. 2: (color online) (a) Contribution of the Mott-insulator 
state to the ground state pmi = K'i/'IMI}] 2 of (|37p as function 
of the ID on-site interaction strength in units of End (d is 
the lattice constant) (b) Logarithmic difference of pmi, cal- 
culated with the groundstate of (|37[) and ()33|) . (c) Logarith- 
mic difference of the adiabatically eliminated photon number 
{ip\a a\ip) with a from (|30p for the two different ground states. 
The parameters are k — l/y2ujit,rj = 2.35u>r and A' c = — k 
(red, solid line), « = Aur,ti — 12.5a;ij and A' c — k (blue, 
dashed line) and k = AuiR,r) = IGujr and A' c = — k/10 (black, 
dashed-dotted line) . In any of the curves, we set V c \ = and 
Uq = — u)r. Here uir is the frequency corresponding to the 
recoil energy, i.e., Er = h 2 k 2 /(2rn) = hcjR. 
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to zero and neglecting terms with respect to the assump- 
tion (1551), we obtain: 



2 



V 



V2A 



Ql 



0(e 3 ) 



(41) 



This is consistent with the assumption g' 2 ~ O (e 2 )- Here 

we defined A = K-iA' c + iU JB. Putting JUJ into (|40b|) 
and neglecting the terms consistent with the order of the 
expansion, such that g\ ~ O (e), it follows that: 



01 



A + rf/A 



(go -Q2) + (e 4 ) 



(42) 



We simplify this expression, q\ rs 77^4 1 (go — £2), which 
is consistent with the above expansion and substitute it 
into (gOal and ([40c) : 



[Hat: Qo] + 2>i£>2 

A -1 (go - 92) + (go - g2)A r 



f?2 — ~TT [H at , Q2] 

in 



lUn 



J N + JB, Q 2 



A 1 (go ~ 02) + (go - Q 2 )A^ 



(43a) 

— 2kq 2 

(43b) 



In order to formulate a master equation for the particle 
variables we have to use the reduced density operator, 
where we trace over the field variables, i.e., g a t — tr(f?) = 
go + £»2 + O (e 4 ). With (|4^a|) and (|43b|) we see that: 



1 



g a t 



[H at , gat] - iU() 



J N + JB, Q 2 



(44) 



As a further approximation, which is also consistent with 
the expansion order of the assumption (|38| , we set (|43b[) 
to zero and neglect [H a t, Q2] and all other terms smaller 
than O (e 3 ) . Then we can express g2 through qq: 



Q2 



2k 



2 r 



A ^o + goA 1 ' 



(45) 



Within this order of magnitude of e we can replace go 
with Qat, leading us finally to the following master equa- 
tion for the reduced density operator of the particle vari- 
ables: 



• 1 

^at — -T [H a t , Qat\ 

in 



(46) 



1 2k 



JqN + JB,(A 1 g a t + g a t A^~ 



Note that this model also contains a damping part, since 
the operator A is not hermitian. Let us investigate this 
damping, by expanding the inverse of A up to first order 
in J, which is consistent with the order of magnitude 
in (|4"6"]) . Hence we replace A -1 and its adjoint in this 
equation by 



A~ 



iA> 



1 - i- 



iAL 



■B 



(47) 



and its adjoint. Since we are restricted on a subspace of 
constant atom number, the Liouville equation reads as 
follows: 



g a t 



1 

ih 



H„ 



(JU 0V ) 2 



nUoy 2 
i 2 + IS! 2 

k 2 ~K. 



2 k 



JB 



B. 



U A'J 2 



A' 2 ' 



B , 0at 



B, Q a 



(48) 



Obviously, the non-dissipative part of this equation 
agrees perfect with our adiabatically eliminated Hamil- 
tonian (|37|) and the structure of the dissipative part is of 
the same Lindblad form as (|3~4")) . Note that an expansion 
of A^ 1 to higher order in J would also provide us the 
correct next-order term of (|37[) plus an extra term in the 
Liouville-equation, which does not correspond to unitary 
time evolution, as described by a Hamiltonian. This con- 
firms the usefulness of the naive elimination method, also 
used in Ref. [H. 



C. Quantum phase transitions in an optical lattice 



In section Till Al we derived two approximate Hamilto- 
nians ([55)) and ([5T|) describing our system of cold atoms 
in an optical lattice. To a large extend they still im- 
plement the well known BH model, but with parame- 
ters controllable via cavity detuning and some additional 
nonlocal interaction terms. Let us now investigate their 
properties in some more detail. One of the key features 
of optical cavities is the feedback mechanism between 
atoms and cavity field. Hence, computations are a sub- 
tle issue, since the matrix elements in the BH Hamil- 
tonian depend on the field amplitude, which itself de- 
pends on the atomic positions. In principle a rigorous 
treatment would consist of calculating the matrix ele- 
ments ([24]) for every photon Fock state and treating the 
parameters of the BH model as operators. To avoid the 
full complexity of such an approach we will first assume 
only a weak dependence of the Wannier functions on the 
mean cavity photon number (a 'a), which allows us to 
proceed analytically. For any set of operating parameters 
we then calculate the matrix elements in a self-consistent 
way replacing the photon number operator by its average 
in the iteration process. Explicitly this is implemented 
by starting from some initial guess jf, 4 0) , J(°), £7(°) 
in the Hamiltonian (|37p . from which we calculate the 
ground state By use of this state we obtain an 

initial mean photon number (ip(°) [a^ 1 a\ip^ oS> ) , with the 
steady-state field operator ([50)) . Now we can calculate 
the matrix elements Jq X \ E^, jW, again leading 
to a new ground state l - !/^ 1 " 1 ) and a new mean pho- 
ton number (^^[a) a[ip^) . Proceeding iteratively, in 
most cases the fixpoint is reached already after very 
few iterations and the system properties are then cal- 
culated with this self-consistent matrix elements. The 
convergence speed decreases near the resonance for the 
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cavity photon number (cf. Fig. [3]), which occurs for 
A c = UqJqN — k, especially for large U . Introducing 
some damping in the iteration procedure easily resolves 
this issue, though. As we mentioned already before, we 




1 2 

&c/(U J N - 

FIG. 3: (color online). Self-consistent photon number in the 
case of four particles in four wells without on-site interaction. 
Parameters are Uq = — ujr and k = ujr. 



restrict the model on a subspace TLn of a fixed total 
particle number TV in an optical lattice of M sites. A 
basis of Hn consists of the states \N, 0, 0, ... , 0), \N — 
1, 1, 0, . . . ,0), . . . , ]0, 0, . . . ,0, N). Since we are interested 
in the quantum phase transition between the Mott insu- 
lator (MI) and the superfluid (SF) state occurring dur- 
ing the variation of certain external parameters, we in- 
vestigate the contributions of these specific states to the 
ground state of the atomic system. The Mott insulator 
state is a product of Fock states with uniform density 
distribution, i.e., 



| MI) 



(49) 



with n — N/M. In contrast, in a SF state each atom is 
delocalized over all sites. It is given by a superposition 
of Fock states, namely of all possible distributions of the 
atoms in the lattice sites, i.e., 



isp) = £ 
fci,fc 2 ,...,fcj\ 



VM W Vki\k 2 \---k M l 



\ki,k 2 , ■ 



(50) 

with X)i=i &i = N. Although the density in the super- 
fluid state is also uniform (ni)sF = N/M and therefore 
equal to the Mott insulator state, its properties are fun- 
damentally different. This manifests especially in the 
spectra and angle dependence of scattered light, pro- 
viding for new, non-destructive probing schemes for the 
atomic phases [13, IH| ■ 

Let us now investigate the influence of the cavity on 
position and shape of the well-known "classical" MI-SF- 
transition [1,0, @. To do so, we compare the two cases of 



a pure quantum field, i.e., V c i — in (|37|) , and a classical 
field (j) = 0) provided by V c i for generating the opti- 
cal potential. We choose rj in such a way, that at zero 
on-site interaction, gw — 0, both potentials are equally 



deep. As depicted in Fig. [4] the influence of the cavity 
strongly depends on the detuning A c . Two contributions 
arise from the quantum nature of the potential. On the 
one hand the potential depth and therefore the matrix 
elements depend on the atomic state. For a classical po- 
tential this is clearly not the case. On the other hand 
the cavity mediates long-range interactions via the field, 
which corresponds to the B 2 -term in (|5T|) . If a poten- 
tial depth near the phase transition point for the quan- 
tum case is associated with some certain average pho- 
ton number n, then n ± 1 are associated with different 
atomic phases. This means that the ground state of the 
quantized cavity field contains contributions of different 
atomic states, each of them correlated with the corre- 
sponding photon number. In this sense photon num- 
ber fluctuations drive particle fluctuations. Depending 
on parameters the former or the latter effect contributes 
more. In Fig.[¥]this is shown for four atoms in four wells, 
where we calculated the occupation probability for the 
Mott insulator pmi = | | 2 and the superfluid state 

p m = | (-01 SF) | 2 for the ground state \ip) of ([37]) as a 
function of the dimensionless one-dimensional on-site in- 
teraction strength giu / '(cLEr) for a purely classical and 
a purely quantum case. For A c — UqJqN = k, photon 
number fluctuations enhance particle fluctuations, shift- 
ing the superfluid to Mott insulator transition to higher 
values of the on-site interaction [Fig. Ufa)]. However, if 
we choose A c — UqJqN = —k, the influence of the atomic 
state on the potential depth exceeds the cavity-mediated 
long-range interactions, strongly shifting the transition 
to lower values of gm [Fig. Hfb)]. Note, that for this be- 
havior, the cavity loss rate must be - although within the 
bad cavity limit - small enough. For larger k the quan- 
tum effects disappear and the ground states for classical 
and quantum potential coincide. 

To correctly address the long-range interactions, cor- 
responding to the B 2 term in ([37)1 . we calculate the 
contribution of the Mott insulator state to the ground 
state of this adiabatic Hamiltonian including and omit- 
ting the B 2 part, respectively. Although, in the situation 
of Fig. HJb) the net effect enhances the phase transition, 
the cavity mediates long-range coherence via B 2 , which 
can be seen by enlarged particle number fluctuations as 
shown in Fig. [5[ Although the effect is not too strong as 
it depends on J 2 is has infinite range and will get more 
important for large particle numbers. 

Finally, we exhibit the transition from a cavity field 
with quantum properties towards a classical optical lat- 
tice. This relies on the assumption that a very bad cav- 
ity should be almost like no cavity and increasing k, 
but keeping the potential depth constant, approaches the 
classical limit. Hence, the effects of the quantum nature 
and feedback of lattice potential should disappear and 
the ground states for classical and quantum potential co- 
incide. The adiabatic eliminated Hamiltonian then has to 
approach the classical Bose-Hubbard Hamiltonian. This 
is shown in Fig. [5] for a system of four atoms in four 
wells, where we simultaneously increase k and rj, keeping 
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FIG. 4: (color online) Cavity influence of the Mott insula- 
tor to supcrfluid transition by means of a comparison of the 
occupation probabilities pmi and psf for a purely quantum 
field, i.e., V c i = 0, and a purely classical field, i.e., rj = 0, 
as a function of the dimensionless one-dimensional on-site in- 
teraction strength gio/(dER). We choose rj such that both 
potentials are of equivalent depth, V = 5.5Er, for zero on- 
site interaction (gio = 0). The quantum and classical case 
is depicted with solid and dashed lines, respectively. In (a) 
we set (Uo, K, rj) = (—1, l/y/2, \/5.5)ujr and A c — UqJqN — k. 
(b) The same as (a) but with A c — UqJqN — —k. 



UoT] 2 /k 2 = —6E R fixed. For every n we calculated the 
value of the on-site interaction giY> , where the contribu- 
tions of the Mott state and the superfluid state to the 
ground state of (J37J) are equal, i.e., |(^|MI)| = |{^|SF)|. 
This is compared with the corresponding value of the 
interaction strength at the same intersection point of 
a purely classical Bose-Hubbard model with a potential 
depth of V c i = —6E R . We see that the transition occurs 
already at a cavity linewidth of only an order of magni- 
tude larger than the recoil frequency, where the deviation 
is small already. Thus one needs quite good resonators 
to see the quantum shift in the phase transition. 
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FIG. 5: (color online) Influence of the long-range interaction 
on the Mott insulator to superfluid phase transition, mediated 
via the B 2 term in (|37[) . The solid line shows the probabil- 
ity for the Mott insulator state as a function of dimension- 
less one-dimensional on-site interaction strength gm/(dEu) 
for a a purely quantum field, i.e., V c i = 0. The dashed line 
corresponds to the probability for the same Hamiltonian, ne- 
glecting the B 2 term. The parameters are the same as in 
Fig. [lb). 
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FIG. 6: (color online) Value of the on-site interaction gw, 
where the contributions of the Mott state and the superfluid 
state to the ground state of (|37[) are equal, i.e., |(t/>|MI}| = 
|(^|SF)|, as a function of k (solid line) for a system of four 
atoms in four wells. Simultaneously we increase 7/, such 
that Uorj 2 /k 2 = —6Er is fixed. Obviously, the correspond- 
ing value at the same intersection point of a purely classi- 
cal Bose-Hubbard Hamiltonian with V c i ~ —6Er is constant 
(dashed line). Parameters are Uo = — lor, k = 4ujr and 
Ac - U N = -k. 



D. Comparison with the full dynamics of the 
master equation 

Using the approximate adiabatic model with elimi- 
nated field we have found important changes in the 
physics so far. Even stronger effects are to be expected 
in the limit of less and less cavity damping and stronger 
atom field coupling. Let us now investigate some first 
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signs of this and test the range of validity of the above 
model in this limit. To do so we have to resort to nu- 
merics and compare solutions of the full master equa- 
tion ([r5]) with the ground states of the adiabatically elim- 
inated Hamiltonian (|37|) . Obviously solving the full mas- 
ter equation is a numerically demanding task. Never- 
theless, by constraining to few atoms in few wells we are 
able to solve the equations and reveal the essential physi- 
cal mechanisms. The limit of the band model description 
is of course reached for atoms coupled strongly to a cavity 
field with only very few photons and no additional clas- 
sical potential V c i present. Here very strong changes in 
the tunneling amplitudes occur whenever a photon leaks 
out of the cavity and reduces the momentary potential 
depth. This leads to strongly enhanced particle hopping. 
For instance, one can think of the situations "one pho- 
ton present" and "no photon present" , where the atoms 
can freely move within the cavity in the absence of an 
external trap. On the other hand one extra photon can 
almost block hopping. Note that in this case the ground 
state atomic configuration can be close to superfluid for 
a low photon number and close to an insulator state for 
a higher photon number. As our matrix elements depend 
only on the mean photon number (a* a), these differences 
cannot be taken into account in an adiabatic model. 

We can explicitly show this behavior by reducing the 
coupling strength Uq, but keeping the average potential 
depth fixed (equal matrix elements) , by means of a higher 
average cavity photon number, which leads to strongly 
reduced photon number fluctuations. The most simple 
situation to discuss this issue is one atom loaded in a lat- 
tice consisting of only two wells. Here, \l) (|r)) means the 
left (right) of the two wells, with a potential minimum 
at s = (x — 7r). The hopping operator B then de- 
scribes tunneling from the left well to the right well and 
vice versa. In Fig. [7] we show this tunneling behavior 
by plotting the mean position of the single atom (kx(t)}. 
The atomic ground state of this system is the symmet- 
ric state \ipo) = (10 + \ r )) /v2 having a mean position of 
{kx)jjj — 7r/2. Decreasing Uq, increasing r\ and adjust- 
ing A c , yields different mean photon numbers (a^a), but 
equal average lattice potential depth V = Uo{a^a). (We 
do not consider an additional classical potential here.) If 
only few photons are present, we observe large fluctua- 
tions of the field and the system damps fast to the ground 
state. As the photon number increases, the potential ap- 
proximates a classical potential as expected, where there 
is no dephasing. The (nearly) equal oscillation frequen- 
cies show that the matrix elements coincide for the dif- 
ferent photon numbers. This is an interesting feature 
corresponding to the quantum nature of the potential. 
In contrast to the Bose-Hubbard model for a classical 
optical lattice, lattice depth and interaction strength are 
not the only important system parameters. Quantum 
fluctuations of the potential are an additional source of 
atomic fluctuations, playing an essential role in the evo- 
lution of the system. Obviously, if only an external po- 
tential is present and the atom is no longer coupled to 



the cavity field (Uq = 0), the system is undamped, due 
to the lack of the only dissipation channel present, cav- 
ity loss. In this case the Hamiltonian ([2"5|) reduces to 
H = (E + JV c i)B + U/2C, and the atom, initially not in 
the symmetric state, oscillates between the left and right 
well. Note that a more rigorous treatment of operator- 
valued matrix elements - as described in the previous 
section - would be capable of describing this behavior 
correctly. Alternatively for few atoms Monte-Carlo wave 
function simulations of the full Hamiltonian could be per- 
formed, allowing for processes, where the particle leaves 
the lowest band [55| . 




FIG. 7: (color online) Mean position (kx(t)) of a single atom 
in two wells. We adjusted Uo,r/,A c in such a way, that 
the mean number of cavity photons increases, but the lat- 
tice depth stays nearly constant: V = KUo{a) a) = —8Er. 
Starting with (Uo,rj,n) — (—50,10,25) (in units of ujr) and 
A c = Jof/o, followed by a successive reduction of Uo by a fac- 
tor of 5, together with an increase of 77 by a factor of y/E and a 
proper adjustment of A c , this leads to mean photon numbers 
of 0.16 (solid line), 0.8 (dashed line), and 4 (dashed-dotted 
line). Initially, the atom is in the right well. 



Obviously, this enhancement of atom fluctuations for 
low photon numbers also affects the dynamics of several 
atoms. We demonstrate this for the case of two atoms 
in two wells. We assume strong coupling with few cavity 
photons and a strong on-site interaction, which - in prin- 
ciple - inhibits tunneling and drives the system deeply 
into the Mott insulator regime. However, starting from a 
state slightly perturbed from the ground state of the adi- 
abatically eliminated Hamiltonian (|57|) . the system does 
not evolve towards this Mott-like ground state but to- 
wards some other, drastically different state. Increas- 
ing the photon number, while keeping the lattice depth 
constant, reduces the atom fluctuations and keeps the 
system near its adiabatic ground state. This is shown 
in Fig. [5Ja) , where the probability for the system being 
in the Mott insulator regime pmi = |V'mi(^)| 2 i s plotted. 
Again we observe that, the larger the intracavity photon 
number is, the more the potential approaches a purely 
classical one and the more significant the ground state 
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probabilities of ([37]) are. Hence we see that including 
the photon number fluctuations strongly suppresses the 
Mott insulator state by allowing the particles to hop dur- 
ing photon number fluctuations. This is also a strong re- 
striction for the use of our adiabatic model Hamiltonian, 
where only average photon numbers enter the model pa- 
rameters. 

Clearly, some added external classical potential dimin- 
ishes this problem as it can ensure the existence of a 
bound state, independent of the number of cavity pho- 
tons, giving an upper limit to the hopping rate. This is 
demonstrated in Fig. [5Jb) , where a classical potential of 
V c i = —5Er is added. Here for (a^a) = 1.44 the devi- 
ations from the adiabatic ground state are of the same 
order as for V c i = for (at a) = 4.8 [Fig. |5{a)]. Never- 
theless, for not too leaky cavities (k is in an intermediate 
regime) , the regime of validity of the adiabatically elimi- 
nated Hamiltonian ([3"T|) is limited to case where either a 
large purely classical potential or a large photon number 
is given. 

Finally, we investigate the other limit of validity, where 
a rather large external classical potential, but only a very 
low photon number is given, i.e., a weakly driven cavity. 
Here the ground state properties of our model resemble 
to a very high degree those of the ordinary Bose-Hubbard 
model. As mentioned above, an atomic ensemble inter- 
acting with a purely classical potential, has no channels 
of dissipation in the absence of spontaneous emission. So 
unless we prepare the system in its groundstate, it will 
show undamped oscillation. In strong contrast the cou- 
pling of the atoms to an even small intracavity field with 
a very low photon number opens a dissipation channel. 
Although the enhancement of atom number fluctuations 
due to fluctuation induced tunneling is small, this damp- 
ing still can drive the system into a steady state, very 
closely to the adiabatic ground state of ([3~7| . This is 
shown in Fig. [9] for the case of two atoms in two wells. 
Here we prepare, for different values of on-site interac- 
tion, the atoms in a state perturbed from the ground 
state of (|37[) with initially no photon in the cavity and 
a given value of the classical potential V c i = —WEr. 
For giD — 0, the ground state is the superfluid state, 
so Fig.[9](a) is the generalization of Fig. [7] to two atoms. 
Although the photon number is only (a'a) = 1.3 x 10~ 4 , 
the system is driven into its ground state. For increasing 
interaction strength, the Mott insulator state becomes 
more and more favored. Still, the interaction with the 
tiny intracavity field enables damping of the atomic evo- 
lution towards a steady state, very close to the adiabatic 
ground state. 

This leads to the conclusion, that, although the cav- 
ity field may not lead to significant modifications of the 
ground state of the system, the cavity is a useful tool for 
faster preparing a system of atoms in its ground state by 
opening a dissipation channel, so that it decays towards 
an eigenstate of the adiabatically eliminated Hamilto- 
nian ([3"7| . 
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FIG. 8: (color online) Probability of the Mott insulator state 
(*) | 2 f° r two atoms in two wells. Parameters and proce- 
dure as in Fig. [7] but due the second atom the photon num- 
bers are increased. The on-site interaction is U — 0.32Er. (a) 
Vci = 0. The curves correspond to a mean photon number of 
0.19 (solid line), 0.97 (dashed line), 4.8 (dashed-dotted line) 
and 24.2 (dotted line), (b) V c i = —5Er and corresponding 
photon numbers of 0.05 (solid line), 0.28 (dashed line), 1.44 
(dashed-dotted line). 



IV. ATOM PUMPING 



Let us now return to our starting Hamiltonian ([26)) and 
consider a second generic model, where the pump laser 
is not injected through the cavity mirrors, but directly 
illuminating the atoms. This rather small change has a 
drastic influence on the physical behavior of this system. 
In the case of cavity pumping, all atoms are simultane- 
ously coupled to the same mode. In this way the cavity 
field depends on the atomic distribution and long range 
order interactions are mediated via the cavity field, influ- 
encing the Mott-insulator to superfluid phase transition. 
In the new geometry, only the directly excited atoms co- 
herently scatter photons in the cavity mode. Due to the 
position-dependent coupling, the scattered field ampli- 
tude and phase for each atom is strongly position depen- 
dent. Atoms located at nodes are not coupled to the field, 
leading to no scattering, whereas atoms at antinodes 
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FIG. 9: (color online) Probability of the Mott insulator state 
|^mi(£)| 2 for two atoms in two wells for different on-site in- 
teraction. In (a) there is no interaction, i.e., U = 0, in (b) 
U = 0.0065E R , in (c) U = 0.0324£ fl and in (d) U = 0.081£ fl . 
Other parameters are (Uo, k, r/, A c ) = (—50, 25, 1, 0) (in units 
of ujr), the classical potential is V c ; = — WEr. The solid 
hne in each subplot shows the corresponding ground state 
probability of (|37p and the number of cavity photons is 
(a f a) = 1.3 x 10" 4 . 



are maximally coupled, leading to maximum scattering. 
Atoms in adjacent wells are separated by half a wave- 
length and scatter with opposite phases, such that their 
contributions to the scattered field interferes destruc- 
tively. Naively one would thus immediately conclude that 
atoms forming a state with a homogeneous density scat- 
ter no field at all so that nothing happens [H], HH . Nev- 
ertheless, fluctuations of the density still can allow for 
some background scattering which should diminish for 
lower temperature. For suitable parameters the corre- 
sponding forces start to reorder the atoms towards a pe- 
riodic pattern of the atoms, where scattering is strongly 
enhanced. This then deepens the optical potential, sta- 
bilizing the pattern in a self -org anizing runaway process, 
semiclassically described in jl8j |. 

At T = quantum fluctuations still can trigger this 
reorganization. To study this effect we assume the co- 
herent pump field to form a broad plane wave propagat- 
ing transversally to the cavity axis (see Fig. [1]) replacing 
cavity pumping. This means that we set rj = and the 
Hamiltonian (|26[) for constant atom number N reads as 
follows: 

H = (E + JV c i)B + h(U J N - A c ) a ] a + —C 

+hU a Ja [ aB + hr) cS (a + a f ) J D. (51) 

Here we introduced the operator D = ^ k {—t) k+1 h k 
describing the difference in atom number between odd 



and even sites. The corresponding Heisenberg equation 
for the cavity field (f2"7)) reads as follows: 



a = ji A c - U (joN + JB*j - /t| a - iij cS J D. (52) 

Consequently the Heisenberg equation for the particle 
operators is: 

b k = (E + JV d - iUoJa^a) (b k+1 + 6 fc _ x ) 

- ir, cS J (a + a f ) {-l) k+1 b k + Un k b k . (53) 

Hence we see that the occupation number difference 
drives the cavity field, which then in turn starts to de- 
phase neighboring atom sites via the first term in the 
second line of Eq. ([53")) . Note that this interesting part 
of the dynamics even survives for deeper lattices when 
J is negligibly small and Jo is of order unity. This will 
discussed in more detail using various approximations be- 
low. 



A. Field-eliminated Hamiltonian 

Adiabatic elimination of the field variables is a bit more 
subtle here as compared to the cavity pump case dis- 
cussed before. The scattering amplitude of light into the 
cavity mode here depends strongly on the atomic posi- 
tions. Hence even small position changes have a large 
influence on the cavity field dynamics. The maximum 
photon number is established when all the atoms are well 
localized at either only odd or only even lattice sites. For 
red atom field detuning this increases the lattice depth 
and forces the atoms into one of two stable patterns, 
where the wells where atoms are located are deeper than 
the empty ones. Hence this changes the translational 
periodicity of the optical lattice from A/2 to A. Such 
bistable behavior was observed by Vuletic and cowork- 
ers and explained in a semiclassical treatment [HI]. 

Let us now turn to a quantum treatment of atoms and 
field. Naive adiabatic elimination encounters a first diffi- 
culty, as the operators B, D do not commute, [B, D] ^ 0. 
Hence this already creates an ordering problem in the for- 
mal steady-state solution of (|52"|) . which gets even more 
difficult when it comes to the replacement of the field 
operators to obtain an effective Hamiltonian (|?T|) . Un- 
fortunately also the second approach used in the case of 
cavity pumping, namely reading off an effective Hamilto- 
nian from the particle operator Heisenberg equation does 
not resolve this problems. Replacing a with the steady- 
state expression in (|53[) leads to a rather complex form, 
so that there is no simple way to find a suitable effective 
Hamiltonian i? a d, with ihb k = [&k,i?ad]- 

Hence we have to resort to the further approximation 
of neglecting the term HUoJa^aB, compared to JV c i. 
This still leaves the most important part of the new 
physics, but reduces the field equation to the form: 
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a = (iA' c — k) a — irj e gjoD. 



(54) 



The steady-state solution of this equation is immedi- 
ately at hand and free of ordering ambiguities of non- 
commuting operators. 



irjcffJo a 

a = -rr, D. 

iA' - K 



(55) 



Also the particle operator equation is much simpler 
within this approximation: 

b k = (E + JVd) (6fe-i + bk+x) 

- zr; cff J (a + at) (-l) k+1 b k + Uh k b k . (56) 

In this form one then can find a well defined effective 
Hamiltonian only containing particle operators. Let us 
thus proceed as in Sec. IIII Al and simply substitute (155)) 
and its adjoint into (|5~lj) . This yields the effective Hamil- 
tonian: 



H ad = (E+ JV cl )B + | 



h J 2 n 2 A' . 
,LJ 'IcB^c jj2 



(57) 



Within first order in J the replacement of the field vari- 
ables in the Liouvillean part of the master equation (|20|) 
in this case does not provide an extra terms to be in- 
cluded in the Hamiltonian. So the effective cavity de- 
cay induced dissipation of the atomic dynamics takes the 
simple and intuitive form: 



£ a d0 : 



% 2 + A'} 



2DgD 



D 2 g 



gD' 



(58) 



Note that Eq. (|57|) with the replacement of the field 
operator by its steady-state expression also leads to the 
same time evolution as induced by (|56|) after symmetriz- 
ing with respect to the field terms. The two approaches 
thus lead to identical predictions, which we will exhibit 
in some more detail in the following. 



B. Self-organization of atoms in an optical lattice 

In this section we investigate the microscopic dynamics 
of self-ordering near zero temperature and compare the 
results of the general model Hamiltonian (|5~lj) and the 
corresponding effective Hamiltonian (|57|) . In order to 
simplify things, we keep the approximation from above 
and neglect JUoa^a in the model, i.e., 



H = (E+ JVd) B - hA' c a) 



hr) cB J D (a 



(59) 

with A' c = A c — UqJqN. Let us point out here, that the 
Hamiltonian in this approximative form is equivalent to 



a Hamiltonian describing ID motion along an optical lat- 
tice transverse to the cavity axis. Such a lattice can e.g. 
be generated by the pump laser itself as it was studied 
in [5y, |57| to investigate the onset of the self-organization 
process[18l. Il9l . l20j at zero temperature. 

Similar to that case, the effective Hamiltonian Eq. |59|) 
for moderate coupling reproduces quite well the results 
of a full Monte-Carlo wavefunction simulation. We have 
checked this for a rather small system of two atoms in 
two wells with periodic boundary conditions. This is the 
minimal system to study self-organization but in general 
sufficient to capture the physics. In this special case the 
operator B simply couples the ordered |11) state to the 
state l/\/2 (|20) + |02)), while the operator D 2 leaves all 
the basis states {|11), 1 20) , |02)} unchanged. It simply 
leads to a relative energy shift. Hence starting from a 
perfectly ordered atomic state (the analog of the Mott 
insulator state) the Hamiltonian part of the time evo- 
lution of the system couples it to the symmetric super- 
position of ordered states. In an adiabatic limit those 
ordered states are correlated with a coherent field ±a in 
the cavity. Thus without damping the evolution would 
simply read: 



|V(*)> = cos (2wi) |11, 0) + ' Sm ^ a;t) (|20, 2a) 

v 2 



|02,-2a)) 
(60) 

where the frequency u> is given by the (E + JV c \)/h. Here 
1 1 1 , 0) is the state with one atom in each well and zero 
photons, whereas |20,a) (|02, — a)) corresponds to the 
state with both atoms in the left (right) well, and the 
cavity field being in a coherent state with amplitude 2a 
(—2a). The factor 2 is due to constructive interference of 
the fields, scattered by the ordered atoms. In the Mott 
state, the scattering fields cancel each other. 

Note that such an entangled superposition of different 
atomic states and fields cannot be reproduced by any 
classical or mean field evolution and requires a genuine 
quantum description. If on-site interaction is added the 
amplitude of this oscillations decreases due to extra rel- 
ative different phase changes of the self-ordered and the 
Mott state. 

Of course we now have to add the effect of dissipation 
via cavity loss. We will see that even single cavity pho- 
ton decay events strongly perturb the system evolution. 
This can be immediately seen by applying the photon 
annihilation operator to the entangled atom-field state, 



\if)(t)') cx a\ip(t)) cx |20, 2a) - |02, -2a). (61) 

This procedure projects out the Mott contributions to 
the state as they are connected to zero photons. Sur- 
prisingly in addition it also blocks further tunneling by 
introducing a minus sign between the two ordered states. 
At this point coherent atomic time evolution stops un- 
til a second photon escapes and re-establishes the plus 
sign. This then allows tunnel coupling back to the Mott 
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insulator state again. In this sense self-ordering is an 
instantaneous projective process here, where the cavity 
acts as measurement apparatus asking a sort of yes/no 
ordering question. 

The fact, that for transverse pumping the adiabatic 
field state associated with the Mott insulator is an in- 
tracavity vacuum decouples this state from further dy- 
namics even in the presence of dissipation. This creates 
numerical difficulties and prohibits an approximation of 
the dissipative dynamics by the adiabatic ground state 
values of (|57[) only. As soon as a photon leaks out of 
the cavity, the contribution of the Mott-insulator state 
is canceled, no matter how large it, corresponding to a 
given on-site interaction, might be. Hence, every initial 
state evolves into a superposition of the ordered states 
and the ground state values of the effective Hamilto- 
nian do not make much sense. Nevertheless, including 
the damping via the effective Liouvillean Q58|) approx- 
imately reveals the complete dynamics. In Fig. [10] we 
show the results of a Monte-Carlo simulation of the dy- 
namics of the Mott and the superfluid contribution, cor- 
responding to (|59p and compare it with a solution of 
the master equation, consisting of the Hamiltonian (|57[) 
and Liouvillean ([55]) . where the field variables are elimi- 
nated. Furthermore, the restriction of the Hilbert space 
to the two states of (jrJTj) and 1 11, 0) , allows for a proof 
of the accuracy of our assumption, concerning the fast 
evolution of the cavity field. We use the coefficients 
c(t), c(t) (calculated with the Monte Carlo simulation) of 
\^(t)) = c(t)|ll, 0)+c(i)(|20, 2a)±|02, -2a)) to construct 
a purely atomic state \<p{t)) = c(t)\ll) + c(t)(|20) ± |02)). 
Then the mean photon number, calculated with the ef- 
fective photon operator (|5"5|) agrees very well with the 
real mean photon number, i.e., 

A ?ff 2 M*)i#V(*)) « m)Wa\m)- (62) 



V. CONCLUSIONS 

Based on an approximative Bose-Hubbard type model 
descriptions, we have shown that quantum characteris- 
tics of light fields generating optical potentials lead to 
shifts in quantum phase transition points and play a de- 
cisive role in the microscopic dynamics of the transition 
process. While many physical aspects can be already 
captured by effective Hamiltonians with rescaled param- 
eters, cavity mediated long-range interactions also play 
an important role and add a new nonlocal element to op- 
tical lattices dynamics for atoms. In that context even 
small modifications in the setup, from cavity pump to 
transverse pump, have a drastic influence on the behav- 
ior of the system on a microscopic level. We have seen 
that the Bose-Hubbard Hamiltonian for the former sys- 
tem can, in a certain parameter regime, be significantly 
simplified by adiabatically eliminating the field variables. 
Although the cavity has influence on its shape, the Mott 
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FIG. 10: (color online) (a) Contribution of the Mott insu- 
lator state pmi in a system of two atoms in two wells. The 
solid line shows the results of a Monte-Carlo simulation, corre- 
sponding to (I59|) with dissipation via cavity loss. The dashed 
line depicts the solution of a master equation with effective 
Hamiltonian (I57p and Liouvillean (|58|l . The constant line, 
shows the ground state value of this contribution of the effec- 
tive Hamiltonian. (b) The same for the contribution of the 
superfluid state psf- (c) shows analogue results for the mean 
photon number. Parameters are V c \ = — WEr, (k, Uo, rfeff) = 
(4, -0.1, 1) lur and A c = U J N + k. 
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insulator to superfluid phase transition occurs similar to 
classical optical lattices. For transverse pumping this is 
not the case. Here, the fields scattered by the atoms in 
the uniform Mott state cancel and completely suppress 
scattering. In parallel new ordered states with maximal 
coupling of pump and cavity field appear and the dy- 
namics favors a superposition of these two ordered states 
correlated with coherent field states with phase difference 
7T. Hence the dynamics generates strong atom field en- 
tanglement and large effective optical nonlinearitics even 
in the limit of linear weak field scattering. 

Of course the various approximations used to derive 
our effective Hamiltonians still leave a lot of room for 
improvements and we could only touch a very small part 
of the physical effects and possibilities contained in these 
model. Fortunately the experimental progress in this 
field is spectacularly fast and several groups now have set 



up optical lattices with cavity fields [3(1 H3, [H, [H, H(| 
and intriguing potential applications of such systems 
were already proposed [EH , so that one can expect a fast 
and exciting further development of this field. 
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